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Motivation for Current Work 

• NASA N+3 Combustor Development Goals 

— Reduce NOx emissions to 80% below ICAO CAEP6 
standards 

— Enhance state-of-the-art for alternate fuels in small core- 
combustors at higher T3 (950K) and higher Operation 
Pressure Ratios (OPR > 50) 

— Leverage N+2 technology achievements of NASA's 

Environmentally Responsible Aircraft (ERA) project (reduce 
NOx emissions to 75% below ICAO CAEP6 standards) 
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Approach for Current Work 

• Lean-Direct Injection (LDI) concepts being studied by OEMs 
and several injector manufacturers 

• Potential to reduce NOx by enhanced mixing, lean burning in primary 
combustion zones near combustor face 

• All primary air comes into primary combustor zone, no dilution air is used 

• Support N+3 Combustor Development with assessment of 
computational models available in the National Combustion 
Code (NCC) for a candidate LDI-3 injector design 

— Reacting flow (comparisons of Effective Area, Combustor 
Temperature, NOx, CO and Unburnt Hydrocarbons with LDI-2 
experimental data) 
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Methods of Current Work 


• Derive a candidate 'small core' injector/combustor 
configuration from N+2 combustor using sets of 
Airblast and Simplex Injectors split into multiple fuel 
stages 

• Use updated physical models in the NCC to predict 
performance and emissions profiles for 

— a candidate LDI-3 geometry configuration at 
'medium-power' conditions 

— comparison of RANS (non-reacting and reacting) 
and TFNS/VLES (non-reacting) 
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LDI-2 Experimental Configurations 


CONFIGURATION 

PILOT (Airblast) 
(OAS/IAS) 

MAIN 1 
(Simplex) 

MAIN 2 

(AirBlast) 

(OAS/IAS) 

MAIN 3 (AirBlast) 
(OAS/IAS) 

SElement Recess Config9 

57CCW / 57CW 

45CW 

45CCW / 45CW 

45CCW / 45CW 

‘Flat Dome’ ConfiglO 

Simplex 55CCW 

45CCW 

45CW / 45CW 

45CW / 45CW 


‘Baseline” or ‘Flat Dome’ - Exit plane of 
the venturi for all thirteen injectors is flush 
with the main combustor dome 
- Single configuration computed with the 
NCC, compared with NASA GRC data 
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LDI-3 Candidate (derived from LDI-2) 
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Typical Geometry for WFST LDI-3 Array 




M1 M2 


M2 


M1 

Ml 

P 

M2 


M2 Ml 


M1 Simplex 
M2 Airblast 
P Simplex 


10M element all-tetrahedral mesh 
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Overview of LDI-3 Injection Elements 

Axial-Bladed-Swirlers with Converging-Diverging Venturi 

Simplex or Airblast Fuel Nozzle 



Various combinations of Simplex and Airblast Elements can be used to form an 
ARRAY of injectors to achieve operability, efficiency and emissions targets 
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Spray Initialization for Simplex Injectors 



Simplex: 0=60'' d0=lO" a=0" (3=0" 
60" ‘hollow’ cone of 10" thickness 



- 32 ‘streams’ per injector 

1 0 particle groups per stream 
Injection location and velocity can be varied 
stochastically as computations proceed 


dn 

n 


4.21x10^ 


d 


n3.5 




, 0.4 


-16.98 


V‘^32/ 


dd 


d 


32 


d 


32 


sauter mean dia = d ^2 < 10fim ; dd=0.2/nm 
number of droplet groups = 8 

dn = number of droplets in the size range d and d + dd 


Y. El Banhawy and J.H.Whitelaw, Calculation of the Flow Properties of a Confined Kerosene-Spray Flame, 
AIAA J., Vol. 18, pp. 1503-1510, 1980 
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Spray Initialization for Airblast Injectors 



Airblast : 0=10° d0=O° a=0° |3=0° 
10° ‘solid’ cone 


- Airblast ‘film injection’ modeled with 8 or 16 
discrete holes for each injector 

- 8 ‘streams’ per discrete hole 

1 0 particle groups per stream 
Injection location and velocity can be varied 
stochastically as computations proceed 


M.S. Raju, LSPRAY-IV: A Lagrangian Spray Module, NASA CR-20 12-2 17294 
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Reduced Mechanism for Jet-A Surrogate 

• 14-species, 18-step finite-rate chemistry model (Ajmani et al AIAA2010-1515) 

• Jet-A surrogate chemistry, mixture of decane (73%), benzene(18%), hexane(9%) 

• Adiabatic flame temperature, flame-speed, ignition-delay matched with shock-tube data 

• allows for in-situ, coupled, computation of emissions (NOx, CO) 


No. 

Reaction 

A 

n 

E 

1 

cllh21 -P o2 => llch + lOh -P o2 

1.00E-P12 

0.00 

3.10E-P04 


forward /cllh21 0.8/ 





forward / o2 0.8/ 




2 

ch -P o2 => CO -P oh 

2.00E-P15 

0.00 

3.00E-P03 

3 

ch -P 0 => CO -P h 

3.00E-P12 

1.00 

O.OOE-POO 

4 

h2 -P o2 < = > h2o -P 0 

3.98E-P11 

1.00 

4.80E-P04 

5 

h2 -P 0 < = > h -P oh 

3.00E-P14 

0.00 

6.00E-P03 

6 

h -P o2 < = > 0 -P oh 

4.00E-P14 

0.00 

1.80E-P04 

7 

h2o -P o2 < = > 2o -P h2o 

3.17E-P12 

2.00 

1.12E-P05 

8 

CO -P oh < = > co2 -P h 

5.51E-P07 

1.27 

-7.58E-P02 

9 

CO -P h2o < = > co2 -P h2 

5.50E-P04 

1.28 

-1.00E-P03 

10 

CO -P h2 -P o2 < = > co2 -P h2o 

1.60E-P14 

1.60 

1.80E-P04 

11 

n -P no < = > n2 -P 0 

3.00E-P12 

0.30 

O.OOE-POO 

12 

n -P o2 < = > no -P 0 

6.40E-P09 

1.00 

3.17E-P03 

13 

n -P oh < = > no -P h 

6.30E-P11 

0.50 

O.OOE-POO 

14 

n -P n -P m < = > n2 -p m 

2.80E-P17 

-0.75 

O.OOE-POO 

15 

h -P n2o < = > n2 -P oh 

3.50E-P14 

0.00 

7.55E-P02 

16 

n2 -P o2 -P 0 < = > n2o -P o2 

1.00E-P15 

0.00 

3.02E-P02 

17 

n2o -P 0 < = > 2no 

1.50E-P15 

0.00 

3.90E-P04 

18 

n2o -P m < = > n2 -P 0 -P m 

1.16E-P15 

0.00 

3.32E-P04 



Use A (pre-exponential factor), n (temperature exponent) and E (activation energy, cal/mol) to compute the Arrhenius rate 
coefficient, k = A(T/Tq)" for a given temperature, T (K). R = universal gas constant, Tq (K) is a reference temperature. 
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Typical Staged NCC Computational Procedure 

• Non-reacting CFD solution 

- k-8, variable Cp, non-linear, cubic, with pressure-gradient effect and wall functions 

- Inlet BC: specify mass-flow rate, static temperature ; Exit BC: specify static-pressure ; All 
walls treated as adiabatic walls 

• Reacting Flow Solution 

— Use ignition sources in region downstream of venturi-exit - ignition sources 
turned off once temperature in any element in mesh reaches 1600K 

- Use spray parameters (SMD, velocity) provided by Woodward FST 

- Lagrangian spray model (spray angle, 0 = 60, spray thickness angle d0 = 10, 32 
streams, 8 droplet groups, stochastic model, no secondary breakup) 

- Finite-rate chemistry models - 14species, 18 steps (direct NO computation) 

• Account for prompt, thermal, N20 NOx pathways 


ISABE Paper 2015-20245 


12 



Validation of NCC RANS for LDI-2 Configurations 
Predicted vs Experimental EINOx Data 


NCC vs Experiment (EINOx) 



NCC Prediction (EINOx) 


13 


RANS vs TFNS/VLES 
Non-Reacting Flow 
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Axial Velocity (m/s) Contours, RANS vs TFNS 

Main 2 AB (Top), Main 1 PA (Bot) (Y=0.034) Main 1 PA (Top), Main 2 AB (Bot) (Z=0.034) 

z RANS 
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Axial Velocity (m/s) Contours, RANS vs TFNS 

+10mm RANS +20mm +114mm 

M2 Ml 


Ml M2 

IM2 p M1 

Ml M2 

c mm -ML. m 

Vax: -20 0 20 40 60 dO 100 120 Vax: -20 0 20 40 60 80 100 120 Vax; -20 0 20 40 GO 80 100 120 

TFNS 

M2 Ml 

Ml M2 

M2 P Ml 

m M2 

d m r. a d m 

Vax: -20 0 20 40 60 80 100 120 Vax: -20 0 20 40 60 80 100 120 Vax: -20 0 20 40 60 80 100 120 
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RANS Reacting Flow CFD 
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Axial Velocity (m/s) Contours, RANS Reacting 

Transverse Cross Sections 


Pilot PA.5.5.CW. (Y=0) RANS Main 1 PA (Top), Main 2 AB (Bot) (Y=0.014) 





M1 


P 



Vax; -20 0 20 40 60 80 100 120 

Main 2 AB (Top), Main 1 PA (Bot) (Y=0.034) z 


M2 


M1 


Vax; -20 0 20 40 60 80 100 120 


Vax: -20 0 20 40 60 80 100 120 



Vax: -20 0 20 40 60 80 100 120 


AC, 

AC, (NCC) 

CFD %Error 

(in^) 

(in^) 


1.40 

1.39 

-0.07% 
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Temperature (K) Contours, RANS Reacting 

Transverse Cross Sections 


Pilot PA5.5.C.W. (Y=0) 


RANS Main 1 PA (Top), Main 2 AB (Bot) (Y=0.014) 



T: 1000 1200 1400 1600 1800 2000 2200 2400 


Main 2 AB (Top), Main 1 PA (Bot) (Y=0.034) 


M2 

5 






/ ■ 


M1 


T: 1000 1200 1400 1600 1800 2000 2200 2400 



T: 1000 1200 1400 1600 1800 2000 2200 2400 


T4 (CEA) 

T4 (NCC) 

CFD % Error 

(K) 

(K) 

(K) 

1723 

1710 

-0.8% 
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NO mass-fraction Contours, RANS Reacting 

Transverse Cross Sections 


Pilot PA.5.5.C.W. (Y=0) 


RANS 











[ ■ 

NO: 5E-05 0.0002 0.00035 0.0005 0.00065 


Main 2 AB (Top), Main 1 PA (Bot) (Y=0.034) 


M2 




M1 




NO: SE-05 0.0002 0.00035 0.0005 0.00065 


Main 1 PA (Top), Main 2 AB (Bot) (Y=0.014) 



NO: 5E-05 0.0002 0.00035 0.0005 0.00065 


~ 1 

NO: 5E-05 0.0002 0.00035 0.0005 0.00065 
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Typical Convergence History 
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T4(CEA)=1651K 
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Medium Power Case 

Mass-weighted average of 
T4 and EINOx across exit 
plane of CFD Domain 
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Axial Velocity (m/s), Temperature (K) RANS Reacting 





Axial Cross Sections 



M2 Ml 

+10mm 

+20mm 


+114mm 

M1 


M2 




M2 

P 

Ml 





Ml M2 





Cz 


■ 


Q 

■ 

Vax: -20 

0 20 40 60 

ao 100 120 

Vax: -20 0 20 40 60 80 100 120 

Vax; -20 0 

20 40 60 80 100 120 


M2 Ml 





M1 


M2 




M2 

P 

Ml 





Ml M2 





c 




L 


T: 1000 1300 1600 1900 2200 2900 

T: 1000 1300 1600 1900 2200 2500 

T: 1000 1300 1600 1900 2200 2500 


ISABE Paper 2015-20245 


23 


NO mass-fraction, UHC mass-fraction RANS 

Axial Cross Sections 


Reacting 


+10mm +20mm +114mm 

M2 Ml 

Ml M2 

M2 p Ml 

Ml M2 


NO: 5E-0S 0.00025 0.00045 0.00065 


NO: 5E-05 0.00025 0.00045 0.00065 


NO: 



SE-QS 0.00025 0.00045 0.00065 


M2 Ml 

M1 M2 
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M1 M2 

: n : n L □ 

C11H21: 0.01 0.03 0.05 0.07 0.09 0.11 C11H21: 0.01 0.03 0.05 0.07 0.09 0.11 C11H21: 0.01 0.03 0.05 0.07 0.09 0.11 
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RANS Reacting Flow Summary 
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LDI-3 'Candidate' - Medium Power Conditions 
Summary EINOx, Exit Temperature, T 4 (K) 



T3 

FAR 

EINOx (Expt) 

EINOx (NCC) 

T4 (Expt) 

T4 (NCC) 


(K) 


(g /kg Fuel) 

(g /kg Fuel) 

(K) 

(K) 

N+3 

950 

0.0234 

6.27 

6.6 

1710 

1723 

N+2 

810 

0.0261 

3.8 

5.4 

1651 

1698 


• Average Exit temperature (T4) predicted by NCC matches experimental data 
to within 1% for candidate N+3 medium-power conditions. 

• EINOx prediction by NCC (6.6) is within 5% error of extrapolated data (6.27) 
from correlation equation of Tacina et. al 

- [Tacina 2014] Tacina, K.M., Chang, C., He, Z.J., Mongia, H.C., Dam, B., and Lee, P., "A Second Generation Swirl-Venturi 
Lean Direct Injection Combustion Concept," AIAA Paper 2014-3434, AIAA Propulsion and Energy Conference, Cleveland, 
OH, July 2014. 

• EINOx extrapolation for N+3 is from Glenn Research Center flame-tube data 
for N+2 configuration 

EINOx = P 3 °-^ (Dp/p) °-^ 

01=0.0081, bj=0.29, 02=0.35, b2=7.15, 03=0.369, b3=7.37 

fi f 2 f 3 are the equivalence ratios for the P, Ml, M2/M3 stages. 
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Lessons Learned and Future Work 

RANS solutions may be very useful as a first-cut to narrow the design 
matrix at medium-power conditions evaluated here; a stepping-stone to 
time-accurate TFNS/VLES computations. 

TFNS Reacting Flow Is prohibitively expensive for design iterations in 
preliminary design phase, particularly with in-situ emissions computations 

NCC RANS predicts EINOx values to within 5% of extrapolated data for 
medium-power N+3 conditions. Additional CFD predictions of N+3 

performance needed for low- and high-power conditions . 

Future work: Evaluate N+3 configuration with reduced-kinetics 
mechanism optimized for emissions 

• AIAA-2014-3662. A Reduced Mechanism for Combustion of Jet-A in LDI Combustor 
CFD Calculations Kumud AJmani; Krishna Kundu; Shaye J. Yungster 


ISABE Paper 2015-20245 



Acknowledgements 

• Subsonic Fixed Wing (SFW) Project at NASA 
GRC 

• Environmentally Responsible Aircraft (ERA) 
Project at NASA GRC 


ISABE Paper 2015-20245 



Backup Slides 


ISABE Paper 2015-20245 


29 



Summary of National Combustion Code (NCC) 

• Two-dimensional axi-symmetric or full three-dimensional computations 

• Most unstructured and structured mesh element types can be computed 

• Finite-Volume solutions of Time-dependent, Navier-Stokes equations 

• Dual time-stepping for 2nd order time-accuracy with 4-stage Runge-Kutta scheme 

• 2-equation, k-8 turbulence models (non-linear, low-Re or wall-functions) (“Generalized 
Wall Function for Complex Turbulent Flows,", T.-H. Shih, LA. Povinelli, K.- H. Chen, N.-S. 
Liu, NASATM 2000-209936.) 

• Lagrangian spray-modeling with primary/secondary breakup and atomization options, 
multi-component fuels (“LSPRAY-IV: A Lagrangian Spray Module”, M. S. Raju, NASA 
CR-2012-217294, Glenn Research Center, Cleveland, OH.) 

• Reduced-kinetics, Finite-rate chemistry models of varying complexity available for various 
fuels 

• Turbulence-chemistry interaction modeled with one of several different approaches 

• RANS time-integration and/or VLES with Time-Filtered Navier-Stokes (TFNS) approach 
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Chemistry Mechanism 


• 14-species, 18-step finite-rate chemistry model (Ajmani et al AIAA2010-1515) 

• Jet-A surrogate chemistry, mixture of decane (73%), benzene(18%), hexane(9%) 

• Adiabatic flame temperature, flame-speed, ignition-delay matched w/experimental 
shock-tube data 

• Iterate with CHEMKIN to find activation energy for fuel breakup step which 
produces a close match between computed ignition delay and experimental ignition 
delay for a particular , P initial temp. (T) 
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Activation Energy vs Temperature 



Paper%15-2024§?° 


1000 

Temperature, K 


Temperature dependent activation 
energy as Jet-A is not a single 
compound 

For CFD code, adjust activation energy 
in breakup step based on initial 
temperature 


1100 


1200 


1300 


31 



